Bose-Einstein condensation and a two-dimensional walk model 
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We introduce a two-dimensional walk model in which a random walker can only move on the first 
quarter of a two-dimensional plane. We calculate the partition function of this walk model using 
a transfer matrix method and show that the model undergoes a phase-transition. Surprisingly the 
partition function of this two-dimensional walk model is exactly equal to that of a driven-diffusive 
system defined on a discrete lattice with periodic boundary conditions in which a phase transition 
occurs from a high-density to a low-density phase. The driven-diffusive system can be mapped to 
a zero-range process where the particles can accumulate in a single lattice site in the low-density 
phase. This is very reminiscent of real-space Bose-Einstein condensation. 

PACS numbers: 05.70Fh, 05.70.Ln,05.50.+q 

Keywords: Walk model, driven-diffusive system, non-equilibrium phase transition 



I. INTRODUCTION 

Over the last couple of years there has been a grow- 
ing interest in studying of the connections between the 
one-dimensional driven-diffusive systems and the two- 
dimensional walk models |l|-|4( . It has been shown that a 
properly defined steady-state normalization factor, called 
the partition function, of some of the one-dimensional 
driven-diffusive systems with open boundaries obtained 
using a matrix product method (reader can see Q for 
a review) is equal to the partition function of a two- 
dimensional walk model obtained using a transfer matrix 
method [J, HQ. 

It is known that the one-dimensional driven-diffusive 
systems exhibit a variety of interesting critical behav- 
iors, such as non-equilibrium phase transition and real- 
space condensation, by changing their microscopic reac- 
tion rates. In Q the authors have shown that the phase 
transitions in the steady-state of non-equilibrium systems 
can be investigated and classified using the Lee- Yang the- 
ory of phase transitions. The asymmetric simple exclu- 
sion process (ASEP) is known as the simplest system 
among the one-dimensional driven-diffusive systems 
In an inhomogeneous version of this process (also known 
as the ASEP in the presence of the impurities or second 
class particles) the particles hop in a preferred direction 
with different hopping rates on a one-dimensional lattice 
with periodic boundary conditions (9141^]. It has been 
shown that there is an exact mapping from this process 
to a zero-range process in which the particles behave 
like bosons and the steady-state is a product measure 
The zero-range process was first introduced into 
the mathematical literature as an example of interact- 
ing Markov processes. By fine-tuning of the microscopic 
hopping rates, a finite fraction of a system's mass can 
accumulate within a microscopic region of the system. 
This can be thought of as a traffic jam. In an equivalent 
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zero-range process, it turns out that the average num- 
ber of the particles at a single lattice site in this phase 
can be of order of the system size. This is in contrast to 
the exclusion process where each lattice site can be only 
occupied by a single particle. The condensation in this 
phase reminds us of real-space Bose-Einstein condensa- 
tion [HG1. 

In recent years, different types of walk models have been 
introduced where some of these models can describe 
physical phenomena such as polymer phase transitions. 
Some examples are given in [13j. As we mentioned ear- 
lier, two-dimensional versions of these models have been 
adopted to describe non-equilibrium phase transition 
in one-dimensional driven-diffusive systems with open 
boundaries. To the best of our knowledge the walk mod- 
els, counterpart of the one-dimensional driven-diffusive 
systems with periodic boundary conditions, have not 
been widely studied. 

In this paper we introduce a two-dimensional walk model 
and calculate its partition function using a transfer ma- 
trix method. We then investigate the critical behaviors 
of this walk model. Finally, we show that the partition 
function of the walk model is exactly equal to that of an 
exactly solvable one-dimensional driven-diffusive system 
with periodic boundary condition. It is known that this 
driven-diffusive system can be mapped to a zero-rage pro- 
cess in which a phase transition into a Bose condensate 
occurs. 



II. THE WALK MODEL 

Consider a walk model in which a random walker starts 
from the origin (0, 0) and takes a finite number of steps 
on = : i, J > are integers} according to two 

forthcoming rules. We assign a weight to each step taken 
by the random walker. All the paths made by the random 
walker are weighted. The weight of a given path will be 
equal to the product of the weights of the consecutive 
steps in that path. The random walker moves according 
to the following rules: 
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FIG. 1: Four different paths of length 5 which end up to 
different heights 0, 1, 2 and 3. 



• For i > j from the lattice site (i, j) to (i + l,j + 1) 
1 •' — 0, 1, 2, • • ■ , oo. This is 



with a weight - where i,j 
to be referred to as an upward step. 



• For i > j from the lattice site (i, j) to (i + 1, 0) with 
a weight p J where i,j = 0, 1, 2, • • • , oo. For j ^ 
(j = 0) this is to be referred to as a downward 
(horizontal) step. 

We will be finally interested in the paths with fixed num- 
ber of steps (fixed length) which contain a certain number 
of downward and horizontal steps (equivalently upward 
steps); therefore, for our later convenience we introduce 
an ad hoc fugacity z and change the second rule as fol- 
lows: 

• For i > j from the lattice site (i, j) to (i + 1,0) with 
a weight zp 3 where i, j = 0, 1,2, • • ■ , oo 

while keeping the first rule unchanged. As can be seen 
the random walker does not take any steps in the negative 
j-direction. Let us assume that the random walker starts 
from the origin (0, 0) and takes N — 1 consecutive steps 
according to the above mentioned rules. The reason that 
we have chosen the number of steps as N — 1 will be clear 
later. After taking these steps the random walker can 
get to the lattice site (N — 1, j) where j = 0, 1, • • • , N — 
1 through different paths. In FIG. [T] we have plotted 
four different paths of length 5 according to the above 
mentioned rules. It is easy to see that there is only one 
path which ends up at the lattice site (TV — 1,N — 1) 
and has the weight l/p w_1 . It can be verified that there 
are 2 N ~ J ~ 2 different ways to get to the lattice site (N — 
for j = 0, 1, • • • , N - 2. The total weight of the 
paths which end up at the lattice site (N — for j = 
0, 1, • • • , N - 2 is equal to z(z + l) N ~ 3 - 2 /p>. 
We define the partition function of the walk model as 
the sum of the unnormalized weights of different paths 
consisting of N — 1 steps. In what follows we show that 
the partition function of the model can be obtained using 



a transfer matrix method. We assume that after i steps 
the state (position) of the random walker is given by the 
vector |j) where j is the height of the random walker 
measured from the horizontal axis and can be an integer 
from to i. These vectors have the following properties: 

\j)k = Sj,k for 3, k = 0, 1, • • • , oo, 
01 f) = f OT j,f = 0, 1, • • • , oo, 

in which T is an infinite-dimensional identity matrix. 
Now we define a transfer matrix T with the following 
property at each step: 



T\j) = V|0) + i|i + l). 



(1) 



One can introduce the following matrix representation 
for the transfer matrix T: 
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(2) 



In fact the transfer matrix T updates the state (position) 
of the random walker. The random walker starts from 
the origin (0,0) which is represented by |0) at the zero- 
th step. After taking N—l steps the random walker is in 
the state \j) where j = 0, 1, • • • , N — 1. As we mentioned 
above, there arc different ways (along different paths) 
to get to the lattice site (N — The total weight 

associated with the lattice site (N — is equal to the 
sum of the weights associated with the paths that end at 
that lattice site. It can be verified that: 



^" 1 |0) = E iii± ^ li> + ^rl^-i>. (3) 
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Finally the partition function of the model, which is equal 
to the sum of the weights associated with the lattice sites 
(N — with j = 0, 1, • • • , N — 1, can be easily calcu- 
lated by multiplying J2fLo(j\ by © from the left: 



Z N (p,z) = X^IT^IO) 

N-2 

- E 

3=0 



z(z + l) N -i- 2 



p' 



N-l 



p 



(4) 



N-lN-i—l , 



»=1 3=0 



JV-1 



in which: 



6!(o-6)! 
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is the usual binomial coefficient. As we mentioned above, 
we are interested in the partition function of the original 
walk model in a special case that after taking N — 1 
consecutive steps, the random walker has taken exactly 
N - M upward steps in which 1 < M < N. 
In order to find the partition function of the model in 
this case, let us have a closer look at the role of the 
fugacity z. The weight associated with a horizontal or 
downward movement is proportional to z; therefore, the 
coefficient of z M_1 in (J4j) is equal to the partition function 
of the walk model Zm,m{p) which consists of exactly iV — 
M upward steps. One can easily use (0} to obtain the 



coefficient of 



,M-1 



in Zn(j>, z). For M/l the result is: 
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and obviously for M = 1 one finds: 
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(5) 
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The partition function ([5]) has a simple explanation in 
terms of the weighted paths in our walk model. It is the 
sum of the weights of the paths of length N — 1 which 
precisely contain N — M upward steps (or equivalently 
M — 1 horizontal and downward steps) . 
As a relevant quantity, one can investigate the mean 
height of the random walker h. The probability of being 
at the height j under the above mentioned conditions is 
given by: 



N,M 



(i) = 



i 



J N,M 



(P) 
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Now the height of the random walker averaged over all 
the steps of each walk and over all the walks is given by: 



N-M 

E^ 



N,M 



d\nZ NiM (p) 
dp 



(8) 



It turns out that in thermodynamic limit N — > oo with 
M = Np, the mean height of the random walker is given 
by: 



N(l -j^) for p < 1 - p, 
p ^rr p for P > 1 - P- 



(9) 



As can be seen, there is a phase transition in the thermo- 
dynamic limit from a phase in which the mean height of 
the random walker is of order N to another phase where 
it is a constant. 

In what follows we show that the partition function of 
the walk model (0 is exactly equal to that of a driven- 
diffusive system with periodic boundary conditions. It is 
known that this driven-diffusive system is equivalent to 
a zero-range process in which a real-space Bose-Einstein 
condensation occurs. 



III. AN EQUIVALENT DRIVEN-DIFFUSIVE 
SYSTEM 

In flll | the author introduces a one-dimensional driven- 
diffusive system of classical particles with hardcore inter- 
actions. The system consists of a particle of type A and 
M — 1 particles of type B moving on a one-dimensional 
lattice of length N with periodic boundary conditions. 
The particle of type A hops from the lattice site i to i + 1 
with rate p provided the target lattice site is empty. The 
particles of type B hop from lattice site i to lattice site 
i + 1 with rate 1 provided that the target lattice site is 
empty. If an empty lattice site is represented by 0, then 
the dynamical rules are simply as follows: 



.4 
B 



A with rate p 
B with rate 1. 



A similar model has also been introduced in [9[ and stud- 
ied in details in [lfj. In this model the particle of type 
A is called an impurity while the particles of type B 
are called the normal particles. The normal particles 
can overtake the impurity. The model can be exactly 
solved using a matrix product method. It turns out that 
the phase diagram of the model consists of four different 
phases. In one of the phases the model presents a shock, 
i.e. a sharp discontinuity between a region of high den- 
sity of normal particles and a region of low density. 
It is shown that the probability distribution of the above 
mentioned system (the one without overtaking) can also 
be obtained using a matrix product method [TlJ. Let 
us label the particle of type A with 1 and the parti- 
cles of type B with 2, 3, ■ ■ • , M — 1 respectively. Accord- 
ing to the matrix product formalism, in the steady-state 
the probability of finding the system in the configuration 
{n} = {77,1,77,2, • • • , Um} in which m empty lattice sites 
lie in front of the ith particle is given by: 



P({n}) 



Zn,m(j>) 



Tr(D' E ni DE" 2 ■ ■ ■ DE nM ) (10) 



in which the operators D', D and E are associated with 
the presence of a particle of type A, a particle of type B 
and an empty lattice site respectively. The denominator 
in (fTU|) is the normalization factor (which is called the 
canonical partition function) and should be calculated 
by considering the conservation of number of empty lat- 
tice sites i.e. Yli=i ni — N — M. The expression (fT0|) 
describes the steady-state of the system provided that 
operators D', D and E satisfy the following quadratic 
algebra [111 ]: 



pD'E = D' 
DE = D. 



(11) 



The partition function Zn,m{p) in (TTOj) has been calcu- 
lated in and it turns out that it is exactly equal to 
that of our walk model given in ([S]) . In the same reference 
the critical behaviors of this driven-diffusive system have 
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been investigated. In the thermodynamic limit L — ¥ oo 
and by defining the density of particles p as M = Np the 
system undergoes a phase transition from a low-density 
phase to a high-density phase. In the low-density phase 
p < 1 — p the number of the empty lattice sites in front 
of the particle of type A is of order N. In [12| the author 
has shown that the above mentioned driven-diffusive sys- 
tem can be mapped to a zero-range process by converting 
the particles into boxes and the empty lattice sites into 
particles. The low-density phase in the driven-diffusive 
system now corresponds to a Bose-Einstein condensate 
phase in the equivalent zero-rage process. In contrast, in 
the high-density phase p > 1 — p the empty lattice sites 
arc uniformly distributed between the particles [TTJ] . It is 
now easy to compare these results with those presented 
in ([5]) for the walk model. It turns out that the mean 
height of the random walker h is equivalent to the num- 
ber of empty lattice sites in front of the particle of A in 
the driven-diffusive system. 

Now we explain (from a mathematical point of view) why 
the partition function of the walk model obtained us- 
ing the transfer matrix method is equal to that of the 
driven-diffusive system obtained using the matrix prod- 
uct method. Let us start with the driven-diffusive sys- 
tem and the matrix representation of its quadratic alge- 
bra. It can be easily verified that the following infinite- 
dimensional matrix representation satisfy the algebra 
(ED: 



2?' = X>><i|, 

2=0 
oo 

£> = 5>'|o><;|, 

2=0 
1 oo 

E = -J2\i + l)(i 



(12) 



in which 



6ij for i,j 



0, 1, • • • , oo. Another ma- 
trix representation for the algebra (fTTj) is given in Q. 
Although the number of particles in the driven-diffusive 
system does not change by the dynamical rules, it is much 
easier to calculate the grand canonical partition function 
of the driven-diffusive system by defining z as the fugacity 
of the particles of type B. Later we will fix this fugacity 
by the density of the particles of type B. For a system 
consisting of a single particle of type A, the unnormal- 
ized weight associated with an arbitrary configuration is 
proportional to Tr(D'X Tl ■ ■ ■ X TL1 ) in which we have 
defined X Ti= o — E and X Ti —\ = D associated with the 
presence of an empty lattice site and a particle of type B 
at the lattice site i respectively. Note that in the canon- 



canonical partition function is given by: 

Z N (p,z) = ^2z M - 1 Tr(D'X Tl ---X TL _ 1 ) 

{n} 

= ^2Tr(D'z^X Tl ---z T ^X TL _ 1 ) (13) 

{Ti} 

= TriD'C 1 *- 1 ) 

in which we have defined C = E + zD. The density of 
the particles of type B is related to their fugacity through 
the following relation: 



z 9 In Zn(v, z) 
lim — . 

n^oo N dz 



(14) 



Using the matrix representation ()12j) , the grand canonical 
partition function (|13j) can be written as follows: 



ical ensemble one has y\ 



L—x 



n 



M — 1. Now the grand 



Z N {p,z) = TrQTlOKilC"- 1 ) 

i=0 

= f^lCMo). (15) 

i=0 

Using (|12p it can be seen that the matrix C is exactly 
equal with the transfer matrix T given in ([%]); there- 
fore, the grand canonical partition function of the driven- 
diffusive system (|13[) is equal to the partition function of 
the walk model (j4| where we define the ad hoc fugacity z 
for the horizontal and downward steps. If we try to im- 
pose particle number conservation and adopt a canonical 
ensemble, then we have to select the coefficient of z M_1 
in (fT5|) which is exactly ([5]) i.e. the partition function of 
the walk model in the case where the number of upward 
steps is equal to N — M . 



IV. SUMMARY AND OUTLOOK 

It seems that the connection between the two- 
dimensional walk models and one-dimensional driven- 
diffusive systems is of a more general validity than a 
couple of examples studied thus far. The ASEP with 
open boundaries is an example which has been studied 
in detail @. In the authors have studied the ASEP 
in the case that its steady-states can be written in terms 
of a superposition of multiple shocks with random walk 
dynamics. It turns out that the partition function of the 
system in this case, is related to that of a walk model 
consisting of multiple Dyck paths. 

In this paper we introduced a two-dimensional walk 
model and calculated its partition function using the 
transfer matrix method. The paths start from the ori- 
gin, but in contrast with the perviously introduced walk 
models, they do not necessarily end on the horizontal 
axis. We showed that this walk model is closely related 
to a one-dimensional driven-diffusive system defined on 
a closed lattice which can be mapped to a zero-range 
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process [TT], LL2J- As it is pointed in [12| a real-space 
Bose-Einstein condensation can occur in this zero-range 
process by fine-tuning of the microscopic hopping rates. 
The zero-range processes have been extensively studied 
in related literatures 14j; therefore, how the zero-range 
processes and the walk models are related still remains 
an open question. 
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